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"Sixth root of unity" and Feynman diagrams: hypergeometric function 
approach point of view 

M.Yu. Kalmykov a B.A. Kniehl a 

a II. Institut fur Theoretische Physik, Universitat Hamburg, 
Luruper Chaussee 149, 22761 Hamburg, Germany 

We briefly discuss the transcendental constants generated through the e expansion of generalized hypergeometric 
functions and their interrelation with the "sixth root of unity." 



1. After the appearance of the proof of the the- 
orem concerning the calculability of 4-loop renor- 
malization group (RG) functions in the frame- 
work of dimensional regularization 1 in arbi- 
trary renormalizable models in terms of £ func- 
tions [2J and the explicit evaluation of some 4- 
and 5-loop RG functions [3], quite intriguing re- 
sults were derived by Broadhurst [4|, who had 
observed the appearance of non-zeta terms in 
the higher-order terms of the e expansion of the 
two-loop massless propagator diagram. Later, 
this observation was clearly explained in the 
framework of the knot approach to Feynman di- 
agrams [5], and all new transcendentals were 
parametrized in terms of multiple zeta values: 

O = E„ 1 >n a >...>n*>oIlJ=l -JT- This result waS 

3 

extended to the case of the 2-loop massless prop- 
agator diagram in [6] and was recently fully an- 
alyzed in [7]. In particular, it was shown that 
the e expansions of 4-loop non-planar massless 
diagrams may generate transcendentals express- 
ible in terms of multiple polylogarithms of the 
"sixth root of unity." All recent results for the 
first few coefficients of the e expansions of mass- 
less propagator diagrams |8] are in full agreement 
with this theorem. Let us recall that the ap- 
pearance in Feynman diagrams of transcenden- 
tal constants related to the "sixth root of unity" 

h yJ2n 1 >n 2 >...>n k >oU 3 = l^-^ , where Pj G 

3 

{0, 1, • • • , 5}, were predicted by Broadhurst in [S] 
in the context of the study of the finite parts of 
massive three-loop bubble diagrams. The lowest- 



weight Broadhurst bases include the follow- 
ing elements: tt 1 m j 2 ln fc 3, Cl 2 (f ) tt* ln J 2 ln fc 3, 
Li 2 (j) tt 1 lir 7 21n" 3, etc., where i,j, k are inte- 
gers (in [9] only constants up to weight 3 were 
presented). Based on the existing partial results 
[TO] , the authors of [UJ pointed out that, in real 
physical diagrams with two massive cuts only, the 
modified set of transcendental constants differ as 
follows: (i) the factor -4= should be present (for 
example, n, In 2, In 3 of [9] should be -4= , In 3 as in 
[UJ); (ii) all elements are generated by values of 
Li a (exp{i6» fe }) and S o ,6(exp{z0fc}), where 



% = f k , 



fc=l,2, 



(1) 



and Li a (z) and S a ,b(z) are the classical and 
Nielsen polylogarithms, respectively; (iii) the ad- 
ditional statement that these constants should 
form an algebra. This construction is in full 
agreement with the all-order e expansions of the 
Feynman diagrams considered in [12] . In [13] , this 
set of constants was extended up to weight 5, and 
it was shown that new transcendentals are gener- 
ated (xs in the notation of [13]) which are not ex- 
pressible in terms of Nielsen polylogarithms. Fur- 
thermore, the appearance of the factor -4= only 
receives a strong explanation in J14I15I16] . This 
factor comes from (i) the evaluation of a special 
type of multiple sums and/or (ii) the structures 
of the coefficients of the all-order e expansions of 
basic hypergeometric functions. In fact, we havqj 



1 The structures of other types of sums and/or hypergeo- 
metric functions were discussed in |16I17I18] . 



reads: 
F 



{l + ai e} p 
l + fe, {l + e l£ } p - 2 



(1 + 2/e) 

2z 



1+1/ 



oo 



lny+JV*fc+i(i/) 



F 



fc=i 

lR+2 /o^^.^.lF-2-R 



{l + a l e} R +\ {2+die} 1 



(l + 2/e) nP _ 2 _ fl (l+c 8 e) 



2,r 



r s=l 



(1+d.e) 



1-2/ 
1+1/ 



lny + ^e fc $ fe+1 (y) 



fc=i 



fc=i 



F 



{l+a t e} K + 3 , {2+d t e} 
-fs, {1+eis} 



P-3-K 



'' L {2 + c l£ } p - 2 - K+L 



(l + d s e) ^ 



2z 



il s=l 



where F is the hypergeometric function pFp-i, 
R, K and L are integers with < R < P — 2, 
< K < P-3, < L < K, the superscripts 
_R and K — L indicate the lengths of the param- 



= , and $k(y),®k(v),$k(y) 



eter lists, y 



and <5>k+i(y) are linear combinations of multiple 
polylogarithms of the square root of unity [TH] of 
weight k, 



**(*/) 



^ Cff t k ln J (y) Li/sx (±y)-Li/s\ (±1) 



(?) 



SJ 



(4) 



Here, c^. are numerical coefficients, s = 
(si, • • • s„) and a = (cti, • • • , a n ) are multi-indices, 
(Tfe = ±1 arc the square roots of unity, j ' + s\ + 
■ ■ ■ + s n = k, and we have used the definition: 



Li/ 



(*) 



mi>---m„,>0 



2 mi fff ••■<" 



The r.h.s. of Eq. (@| corresponds to the Remiddi- 
Vcrmaseren functions [20]. The differential- 
reduction algorithm |15|21|22] allows us to write 
the expansions for arbitrary values of parameters, 
as 



Q{z) P F P - 1 



{Ii+cney 

L + ^ + fe, {Ki + as} 



P-2 



= j2 R i(z)i E <i- ©i+E^wpq- ©] 



(5) 



where Q(z), Ri(z) and Si(z) are polynomials 
and Ii, Ki and L are integers. For two-cut 
diagrams with massive lines, the on-mass-shell 
case corresponds to z = 1/4, so that y = 
exp(±i7r/3), — y = exp(=Fz27r/3) and j-rf = 
T~js, and the e expansions of Eqs. © and §5§ 
have the structures -4= 2feLo e ''^J+ 1 i^k) and 
Efclo £j *i+i+i ( ±6 ^) ' respectively, where $,(z) 
and $j(2) are defined by Eq. (0J, and 9k by 
Eq. (P). 

2. From the all-order e expansions of the hy- 
pergeometric functions constructed in [16] . the 
following set of constants are generated at weight 
k: 

(i) The Remiddi-Vermaseren functions evaluated 
with arguments 9k defined by Eq. (JTJ) , 



£(* 



Li m (±0 fc )-Li m (±1 



V*) 



(?) 



(6) 



This is a subset of the "sixth root of unity" |5], 
which, up to weight 4, agrees with |11) and, at 
weight 5, with [T3]. (ii) The product of Remiddi- 
Vermaseren functions, with arguments 9k defined 
by Eq. (JTJ, multiplied by i/VE, 



71 ^ 

8, J 



Li (l) {±9k)-U (3 ,(±l) 



'(;) 



(?) 



(7) 



Up to weights 4 and 5, these results agree with 
[TTj and [13] , respectively. This explains the mys- 
terious generation of the factor l/-\/3. 

The product of hypergeometric functions can 
be understood as the product of low-loop mas- 
ter integrals (the 2-loop propagator diagram may 
generate the product of two 1-loop self-energies). 
The appearance of diagrams of this type gives rise 
exactly to the product of two basis elements de- 
scribed in [11] , 

For practical applications in the PSLQ anal- 
ysis [23 , it is desirable to have a minimal set 
of constants corresponding to Eqs. ([6]) and (JTJ. 
Their special subclass, Euler-Zagier sums, was 
analyzed in [24]. For the other constants, the 
solution is not unique, since the commonly ac- 
cepted parametrizations of the real and imagi- 



nary parts of the Remiddi-Vermaseren function^ 
with argument z = exp(i0) do not exist. For ex- 
ample, starting from weight 3, the inverse bino- 
■h, can be written (for de- 



mial sums, £) n=1 pjT^ 

tails, see [13 26 27 28 ) in terms of either Dirich- 
let's L series, the derivatives of \& functions, gen- 
eralized log-sin functions [55], Nielsen polyloga- 
rithms [26113] or generalized polylogarithms. Let 
us recall that the classical polylogarithms with 
these arguments produce the Clausen functions 
Clj (9) [33] , while the Nielsen polylogarithms pro- 
duce the generalized log-sine functions Ls^' (0) 
only ]T3 26 29 30 21] . For the parametrizations 
of Remiddi-Vermaseren functions with complex 
unit, hsc a ,b(9) and LsLsc Qi f, !C (9) functions were 
introduced in [14130132] . 

3. The e expansions of hypergeometric functions 
entering the r.h.s. of Eqs.© and ([3]) may be writ- 
ten in terms of multiple inverse binomial sums 
[13114126128] defined as 



^fli ,...,a p 



,b„;c 



S ai ■ ■ ■S ap Ab 1 . . . Aft 



i=i 



(?) 



,(8) 



/ 



where S a and Sb stand for S a (j—1) and S&(2j — 1), 

respectively, Sk(j) — J2i=i l~ k are the harmonic 
sums and A Q denotes the following linear combi- 
nations of S a : 

exp(JTe k S k )=l + eSi + e 2 [S 2 + S 2 ] 

+e 3 [s^ + 3SiS 2 + 2S 3 

+e 4 \S 1 4 + 6S 2 S 2 + 3S 2 + 8S 1 S 3 + 6S 4 \ +0(e 5 ) 

oo 

i=i 

Here, all rational factors coming from the series 
expansion of the exponent are equal to 1 . For the 
hypergeometric functions of Eq. @ , we have 



E 






(?b- ^s 



I I \1)J 



cM K+1) 



(9) 



2 For numerical evaluations of Remiddi-Vermaseren func- 
tions, some of the existing programs | 25| may be used. 



where c r are rational numbers, fli- are products 
of the harmonic sums S a and A&, namely 

n^e{S 2 ,S 2 ,SiAi,A 2 }, 

ni 3} G {S 3 ,S 1 S 2) Sf,S 2 A 1 ,SfA 1 ,S 1 A 2 ,A 3 } , 

fi}, g {S 4 ,SiS 3 , S 2 , S 1 S 2 , Si, A4, Si A3, 

S 2 A 2 ,S 2 A 2 ,SfAi,SiS 2 Ai,S 3 Ai} , (10) 

and ujr belong to sets of transcendental con- 
stants, namely 

u^=ir, 4 2) e{Ls 2 (f),7rln3,}, 
4 3) e{C 3 ,^ln 2 3,7rC 2 }, 

4 4) G {C 4 ,Ls 4 (f) ,7rC 3 ^ln 3 3,C 2 4 2) } • (H) 
At weight 5, there are 10 independent terms: 



4 5) 

,(5) 



C 5 ,w 2 



(5) 



Di , J 3 5) 



.(5) 



Ls 5 

,(5) 



: ir( 4 , w£ ' = ttC 3 In 3 , uj y 6 ' = irC, 2 In 3 , 

^ln 4 3, 4 5) =7r[Ls 2 (f)] 2 , 



,( 5 ) 



,( 5 ) 



C 3 Ls 2 (f ) 



.(5) 
; 10 



C2C3 



(12) 



where Ls^- (9) are generalized log-sine functions 
defined as 



i4 fc) (0) = - 



k \v> 



-fc-i 



2sinf 



Ls, (9) 



I*S 0) (*) , 



(13) 



and C 3 , C 4 , C5 and Di are combinations of gen- 
eralized log-sine functions, namely 

C 3 = 3Ls 3 (f)-2Ls 2 (f)ln3, 

C 4 = 2Ls 4 ( 2 f)-3Ls 3 (^)ln3+Ls 2 (f)ln 2 3, 

C 5 = Ls 5 (f)-2Ls 4 (f)ln3 



Di 



+ -Ls 3 (f)ln 2 3-±Ls 2 (f)ln 3 3 

3Lsf (f)-4.Ls«(f) 

32 

+ 2yLs 4 (f)ln3 + 8C 2 Ls 3 (f) . 



(14) 



We wish to mention that the combinations C 3 , 
C 4 and C5 completely coincide with appropriate 
terms of the e expansions of the Gauss hyper- 
geometric functions considered in [12] , and the 



generating function for these combinations is 



E^ +1 = ^e^ls, +1 ^ 

3=0 3=0 J ' 



so that C\ = and C 2 = 2Ls 2 (§) . For the hyper- 
geometric functions of Eq.Q, the e expansions 
have the form 



E 



n 



(R) 



(?) 



3 =i I 7 I ja+T s=1 



= £C s a< 2 + r 



+fl) 



(15) 



where r is a positive integer, C s is a rational num- 
ber and ov belong to sets of transcendental con- 
stants, namely 

<4 2 >=C2, ^ 3 )e{C 3! 7rLs 2 (f)}, 

aWe{^Ls 3 (f),[Ls 2 (f)] 2 ,C4}. (16) 

At weight 5, there are only 6 independent con- 
stants, namely 

(j[ 5) = 7rLS4 (f ) , 4 5) = 7rLs 4 (^) , 

o-3 5) = ttC 2 Ls 2 (f ) , o-f = X5 . 

&5 — ^5 > a 6 — S.2C.3 j 



(17) 



where xs is defined as %5 = E^i prS [13130] . 
For illustration, we present here the analytical re- 
sults for all sums of weight 5: 



04 



,(5) 



155 



,(5) , 8 ,,( 5 ) 



, ,, 9-3 +"54-«<4-3«8 +3-9 



5|\ _ 4 i(5) 233 



(5) , »,,(5) 



SiS 8 \ 9 107 (5) 2437 

-r/ = _ 8 Ci+ ^ W3 + l32- <4 

, 16 ,,(5) , 2 ,,(5) 14 , (5) 1(B) 

"9" 5 9 8 ~~9 9 ~3 10 ' 
S?S 2 \ 5 131 (5) 1843 

"■' 3 216 



12 81 



7tC 4 



98 



,(5) 



14 f 5 i 68 



..(5) 1* (5) DO ( 5 ) 
-U)g +-—UJ0 —-—UJa 



81 ° 27 D 27 8 27 

2 (5) 



S *\ 1«^ j_ 23 n 218 (^ 2837 , 



716 (5 ) 110 (5) 1 (5) 2 (5 ) 

+ yT 5 + ~ cj6 + 3 W7 ~r 8 

1 16 (5) , 4 (5) 

Y 9 3 w ' 
/5l3Al \-!9r+ 195 D 2717 J 5 > 



40093 
864 



7rC4- 



9 



" W 5 



6' 



13 ,, (5) + 94 M (5 U 5 ,/ 5) 

~T 8 ~9 9 3 10 ' 

' gAAi \ _ 83 1493 ( 5 ) 

v j I 48 1 324 ^ 3 

50435 278 (5) 1 (5) 

7TC4 W H CJr 

2592 s 81 5 27 6 

, 17 (5) 74 ,,(5) 1 (5) 
+ 27^ 8 "27^ + 9 Wl ° ' 
"S 3 Ai\ 9 89 (5) 

8 54 3 



1 (5) 



J I 
1855 „ 

H 7rC4- 

432 



62 (; 
9 9 



9 5 9 8 



:s) 



1(5) 

"3 Wl ° 



3470 (5) 482 (5) 



257 3593 (5) 

162 3 



9' 



102433 > 3470 (5) 482 

+ 1 (5) 56 w (5) 404 (5) 14 (5) 
3 7 ~27 8 ~27~ 9 "o" 10 
/ ^ 2 A 2 \ 73 1249 (5) 

v j I 24 162 3 

42319 458 (5) 1 (5) 

- mT^-sT^ + 27 We 

16 (; (5) 136 (5) 4 (5) 
~^ W 8 + ~27~ W 9 + 9 Wl ° ' 

r+ 61 D 94 V 5 > 

16771 , 1402 (5) 383 (5) 

144 ^ 4 27 5 18 6 

, 1 (5) 16 (6) ,124 (5) 4 (5) 

3 7 9 8 ~9~ 9 3 w 



"27 
l S 1 A 3 \ 



/A4 X 
\ J 



5093 



7 422 f Bl 

28C 5 + -^-— W f 

1576 (5) 226 (5) 

w 5 +"^^6 



■7rC4 



36 " S4 ^ 27 ' y 

1 ,,( 5 ) 16 ,,( 5 )-l 136 , ,( 5 ) + 4 ,/ 5 ) riR^ 



Si 

J 4 

Ai 

J 4 
Ss 

9 

n 

n 3 
Si 02 



28 (6) ,19. . 134 
""Si" 1 + 27 C2C3 + ^7 C5 

82 (6) ,46 847 

^sT 1 + 27 C2C3 + ^T C5 



.(5) 



8 58 

+ qC2C3 + -^"C5 



r ( 5 ) 



r ( 5 ) 



27 



H- 



SaAj 



9 J 
4 
"27 C 

243 4 81 3 

81 C2C,3 gl C,5, 

81 3 



2 29 

+ 36C8+27C1 



.(5) 



146 
243 ( 



95 662 

~^-C2C3-" 8 yC5 



r5| 
S1A1 



2 (5) 1 
= 9 CTl "S^ 



13 

18 



C2C3 



47 
18< 



;?•' 



277 
324 ( 



323 
432^3 



.(5) 



Ci — — (T 



(5) 



23 

108 ( 



1291 11 
144 ^~Tq X5 ' 



r (5) 



5?A! 



116 
243 ( 



.(5) 



,(5) 



23 

81 ( 



.(5) 



437 529 3 

162 C2C3+ 162 Cs+ 2 X5 ' 



337 
162 C 



A3 



216 ^2S3 

1015 



.(5) _(5) 
1 CT 2 

6037 
216 



67 

— c 

54 

7 



.(B) 



<5-oX5 , 



.(5) 



. 7 (5) 469 
324 4 2 2 108 3 
5633 19123 49 

l32" C2C3 + ^32" C5+ 16 X5 ' 



.(5) 



SiA 2 



517 
324 C 



.(5) 



11 

T 



.(5) 



143 
108 f 



.(5) 



3059 7049 35 

l32 _C2 ^ 3 + l32 _< ' 5+ 16 X5 



(19) 



where (X) = ^J2n=i^ X 



[X] 



E 



„--i T^rrA" and C5, £>i, wf and af- are de- 
fined by Eqs. (fT^j) . ([Tl f and (|T7|l. 
4. Based on the recently established analyti- 



cal structures of the coefficients of the all-order 
e expansions of hypergeometric functions, indi- 
cated by Eqs. ([3]), ^ and ([5]), we presented in 
Eqs. dnj, (|nj|, (|TBJ| and (|T7|) the set of linearly 
independent transcendental constants generated 
by single-scale massive Feynman diagrams with 
two massive cuts in 4 — 2e dimensions. The main 
difference between these two sets is the factor -4= , 
which was predicted in [TT]. We mention that 
the first set, given in Eqs. (ITU and (IT2"j) , does 
not form an algebra: to r x u; r 7^ Ur • On 



the other hand, the second set, given in Eqs. ([16]) 
and (IT71) . does form an algebra. Changing the 
space-time dimension or considering Feynman di- 
agrams related to hypergeometric functions of a 
few variables or with a more complicated com- 
bination of massive and massless cuts generates 
another set of constants, as may be seen, for ex- 
ample, in 30 32 33 . 
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